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Abstract 

The exclusive reconstruction of D° mesons in the ALICE experiment 
allows to study the QCD energy loss of charm quarks in the decon- 
fined quark-gluon plasma (QGP) medium expected to be produced in 
central nucleus-nucleus collisions at the Large Hadron Collider. 

1 Introduction 

The ALICE experiment |2j at the LHC will study nucleus-nucleus (AA) col- 
lisions at a centre-of-mass energy -y/s NN = 5.5 TeV (for Pb-Pb) per nucleon- 
nucleon (NN) pair in order to investigate the properties of QCD matter 
at energy densities of few hundred times the density of atomic nuclei. In 
these conditions a deconfined state of quarks and gluons is expected to be 
formed jS]. 

Hard partons and heavy quarks, abundantly produced at LHC energies in 
initial hard scattering processes, are sensitive probes of the medium formed 
in the collision as they may lose energy by gluon bremsstrahlung while prop- 
agating through the medium itself. The attenuation (quenching) of leading 
hadrons and jets observed at RHIC j3j is thought to be due to such a mecha- 
nism. The large masses of the charm and beauty quarks make them qualita- 
tively different probes, since, on well-established QCD grounds, in-medium 
energy loss off massive partons is expected to be significantly smaller than 
off 'massless' partons (light quarks and gluons). Therefore, a comparative 
study of the attenuation of massless and massive probes is a promising tool 
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to test the coherence of the interpretation of quenching effects as energy loss 
in a deconfined medium and to further investigate the properties (density) 
of such medium. 

In the first part of this paper, we shortly summarize one of the widely 
used models of parton energy loss and we discuss how we used it in our sim- 
ulation. In the second part, we show that the exclusive reconstruction of 
D° — > K~7T + decays with ALICE allows to carry out the mentioned compar- 
ative quenching studies by measuring the nuclear modification factor of the 
D mesons transverse momentum (p t ) distribution 

dA^AA/dpt /binary NN collision 
diVpp/dpt 

which would be 1 if the AA collision was a mere superposition of independent 
NN collisions without nuclear or medium effects, and the D/ 'charged hadrons 
(D/h) ratio 

R m (Pt) = R^ A (Pt)/R h AA(Pt)- (2) 

2 Parton energy loss and the dead cone effect 
for heavy quarks 

In this work we use the quenching probabilities (or weights) calculated by 
C.A. Salgado and U.A. Wiedemann [H] in the framework of the 'BDMPS' 
formalism [Oj, which we summarize in the following. The energy loss obtained 
with the quenching weights is presented in Section El 

An energetic parton produced in a hard collision radiates a gluon with 
a probability proportional to its path length L in the dense medium. Then 
(Fig.Q left) the radiated gluon undergoes multiple scatterings in the medium, 
in a Brownian-like motion with mean free path A which decreases as the 
density of the medium increases. The number of scatterings of the radiated 
gluon is also proportional to L. Therefore, the average energy loss of the 
parton is proportional to L 2 . 

The scale of the energy loss is set by the 'maximum' energy of the radiated 
gluons, which depends on L and on the properties of the medium: 

uj c = qL 2 /2, (3) 

where q is the transport coefficient of the medium, defined as the average 
transverse momentum squared transferred to the projectile per unit path 
length, q = (g t 2 ) mcdium /A 0. 
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Figure 1: Typical gluon radiation diagram (left) and transport coefficient 
as a function of energy density (right) for different media: cold (marker), 
massless hot pion gas (dotted curve) and ideal QGP (solid curve) [7j. 



In the case of a static medium, the distribution of the energy uj of the 
radiated gluons (for uj <C uj c ) is of the form: 



where Cr is the QCD coupling factor (Casimir factor), equal to 4/3 for 
quark-gluon coupling and to 3 for gluon-gluon coupling. The integral of the 
energy distribution up to uj c estimates the average energy loss of the initial 
parton: 



The average energy loss is therefore: proportional to a s Cr and, thus, larger 
by a factor 9/4 = 2.25 for gluons than for quarks; proportional to the trans- 
port coefficient of the medium; proportional to L 2 ; independent of the parton 
initial energy. The last point is peculiar to the BDMPS model. Other mod- 
els jHl Ej consider an explicit dependence of AE on the initial energy E. 
However, there is always an intrinsic dependence of the radiated energy on 
the initial energy, determined by the fact that the former cannot be larger 
than the latter, AE < E. 

The transport coefficient is proportional to the density of the scattering 
centres and to the typical momentum transfer in the gluon scattering off 
these centres. Figure [T] (right) reports the estimated dependence of q on the 




(4) 




uj—duj oc a s Cr uj c oc a s Cr q L 2 . 



(5) 
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energy density e for different equilibrated media [Jj: for cold nuclear matter 
(marker) the estimate is g co id — 0.05 GeV 2 /fm; for a QGP formed at the 
LHC with e ~ 50-100 GeV/fm 3 , q is expected to be of ~ 5-10 GeV 2 /fm. 

In Ref. P0| Yu.L. Dokshitzer and D.E. Kharzeev argue that for heavy 
quarks, because of their large mass, the radiative energy loss should be lower 
than for light quarks. The predicted consequence of this effect is an enhance- 
ment of the ratio of D mesons to pions (or hadrons in general) at moderately 
large (5-10 GeV/c) transverse momenta, with respect to what observed in 
the absence of energy loss (proton-proton collisions) . 

Heavy quarks with momenta up to 20-30 GeV/c propagate with a velocity 
which is smaller than the velocity of light. As a consequence, gluon radiation 
at angles smaller than the ratio of their mass to their energy O = m/E is 
suppressed by destructive quantum interference. The relatively depopulated 
cone around the heavy quark direction with < Go is indicated as 'dead 
cone' [TT] . 

In Ref. ^01 the dead cone effect is assumed to characterize also in-medium 
gluon radiation and the energy distribution of the radiated gluons for 
heavy quarks, is estimated to be suppressed by the factor: 



dJ 



Heavy / 



'dJ 

dcu 



Light 



1 + ^ 

e 2 




H/L, 



(6) 



where the expression for the characteristic gluon emission angle jTUj G ~ 
(g/cu 3 ) 1 / 4 has been used. The heavy-to-light suppression factor F h /l in (JUJ) 
increases (less suppression) as the heavy quark energy E increases (the mass 
becomes negligible) and it decreases at large u, indicating that the high- 
energy part of the gluon radiation spectrum is drastically suppressed by the 
dead cone effect. 



3 Simulation of energy loss 

The Salgado- Wiedemann (SW) quenching weight is defined as the probabil- 
ity that a hard parton radiates an energy AE due to scattering in spatially 
extended QCD matter. In Ref. [5], the weights are calculated on the basis 
of the BDMPS formalism, keeping into account both the finite in-medium 
path length L and the dynamic expansion of the medium after the nucleus- 
nucleus collision. The input parameters for the calculation are the length L, 
the transport coefficient q and the parton species (light quark or gluon). 
The distribution of the in-medium path length in the plane transverse to 
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Glauber model for Pb-Pb with b < 3.5 fm (do/db « b) 
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Figure 2: Distribution of the path lengths in the transverse plane for partons 
produced in Pb-Pb collisions with b < 3.5 fm (left). Average energy loss as 
a function of the transport coefficient (right). 



the beam line 1 for central Pb-Pb collisions (impact parameter b < 3.5 fm, 
corresponding to the 5% most central collisions) is calculated in the frame- 
work of the Glauber model of the collision geometry [T2j. For a given impact 
parameter, hard parton production points are sampled according to the den- 
sity Pcoi\(x, y) of binary nucleon-nucleon collisions in the transverse plane 
and their azimuthal propagation directions are sampled uniformly. For a 
parton with production point (xo,yo) and azimuthal direction (u x ,u y ), the 
path length is defined as: 



Many sampling iterations are performed varying the impact parameter b 
from 0.25 fm to 3.25 fm in steps of 0.5 fm. The obtained distributions are 
given a weight b, since we verified that d<r hard /d6 oc b for b < 3.5 fm, and 
added together. The result is shown in Fig. El (left). The average length is 
4.5 fm, corresponding to about 70% of the radius of a Pb nucleus and the dis- 
tribution is significantly accumulated towards low values of L because a large 
fraction of the partons are produced in the periphery of the superposition 
region of the two nuclei ('corona' effect). 

For a given value of the transport coefficient q and a given parton species, 
we use the routine provided in Ref. [Sj to get the energy loss probability 
distribution P(AE; L) for the integer values of L up to 15 fm. Then, these 15 

1 Partons produced at central rapidities propagate in the transverse plane. 
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J °° dl I p C oii(xo + lu x ,y + l u y ) 



(7) 
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distributions are weighted according to the path length probability in Fig. [21 
and added together to obtain a global energy loss probability distribution 
P(AE). The energy loss to be used for the quenching simulation can be 
directly sampled from the P(AE) distribution corresponding to the chosen 
q and to the correct parton species. 

The predicted lower energy loss for charm quarks is accounted for by 
multiplying the P(AE) distribution for light quarks with the dead cone 
suppression factor -Fh/l in (jHJ). It was verified that this approximation is 
equivalent to recalculating the SW quenching weights with the gluon energy 
distribution for heavy quarks as given in (0) [El EH)- Since -F h /l depends 
on the heavy quark energy, the product has to be done for each c quark or, 
more conveniently, in bins of p t . Figure 121 (right) reports the average energy 
loss as a function of the transport coefficient for light quarks and for charm 
quarks (m c = 1.2 GeV) with p t = 1-2, 10, 20 GeV/c, as obtained with the 
described dead cone correction (jo t -dependent P(A.E') ® Fh/l product). With 
q = 4 GeV 2 /fm, our estimated transport coefficient for the LHC (see next 
paragraph), the average energy loss for light quarks is (AE) ~ 35 GeV (the 
effective value of (AE) is lower by about a factor 2, due to the constraint 
AE < E). For c quarks of 1-2, 10, 20 GeV/c (AE) is about 2%, 10% and 
20%, respectively, of the average loss for light quarks. 

For the estimation of the transport coefficient q for our simulation, we 
consider that it is reasonable to require for central nucleus-nucleus collisions 
at the LHC a quenching of hard partons at least of the same magnitude as 
that observed at RHIC. We, therefore, derive the nuclear modification factor 
i?AA for charged hadrons produced at the LHC and we choose the transport 
coefficient in order to obtain Raa — 0.2-0.3 in the range p t = 5-10 GeV/c 
(for RHIC results see e.g. Refs. [SUS]). 

The transverse momentum distributions, for p t > 5 GeV/c, of charged 
hadrons are generated by means of the chain: 

1. generation of a parton, quark or gluon, with p t > 5 GeV/c, using 
PYTHIA jT^] proton-proton with yfs = 5.5 TeV and CTEQ 4L parton 
distribution functions; with these parameters, the parton composition 
given by PYTHIA is 78% gluons and 22% quarks; 

2. sampling of an energy loss AE according to P(AE) and calculation of 
the quenched transverse momentum of the parton, p[ — p t — AE (if 
AE > p t , p[ is set to 0); 

3. (independent) fragmentation of the parton to a hadron using the lead- 
ing order Kniehl-Kramer-Potter (KKP) fragmentation functions |17j . 
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Quenched and unquenched p t distributions are obtained including or exclud- 
ing the second step of the chain. Raa is calculated as the ratio of the p t 
distribution with quenching to the p t distribution without quenching. We 
find Raa ^ 0.25-0.3 in 5 < p t < 10 GeV/c for q = 4 GeV 2 /fm. This value is 
reasonable, as it corresponds, using the plot in Fig. [TJ to an energy density 
e ~ 40-50 GeV/fm 3 , which is about a factor 2 lower than the maximum 
energy density expected for central Pb-Pb collisions at the LHC. 

Charm quarks are generated using PYTHIA, tuned in order to reproduce 
the single-inclusive c (and c) p t distribution predicted by the pQCD pro- 
gram HVQMNR [IE] with m c = 1.2 GeV and /^Fact. = yURonorm. = 2m t = 
2 \/ m c + Pt (t ne details on this tuning can be found in Ref. [H]). We use 
the CTEQ 4L parton distribution functions including the nuclear shadow- 
ing effect by means of the EKS98 parameterization [20J and the parton in- 
trinsic transverse momentum broadening as reported in Ref. ^H]- Energy 
loss for charm quarks is simulated following a slightly different procedure 
with respect to that for light quarks and gluons. Since the total number 
of cc pairs per event has to be conserved, in the cases where the sam- 
pled AE is larger than p t , we assume the c quark to be thermalized in the 
medium and we give it a transverse momentum according to the distribution 
dA/dm t oc m t exp(— m t /T). We use T = 300 MeV as the thermalization 
temperature. The other difference with respect to the previous case is that 
we use the standard string model in PYTHIA for the c quark fragmentation. 



4 Charm reconstruction with ALICE 

The transverse momentum distribution of charm mesons produced at central 
rapidity, \y\ < 1, can be directly measured with ALICE from the exclusive re- 
construction of D° — > K~7T + (and charge conjugates). The displaced vertices 
of D° decays (cr = 124 /im) can be identified in the ALICE Inner Tracking 
System, that provides a measurement of the track impact parameters to the 
collision vertex with a resolution better than 50 /im for p t > 1 GeV/c. The 
low value of the magnetic field (0.4 T) and the K/n separation in the AL- 
ICE Time of Flight allow to extend the measurement of the D° production 
cross section down to almost transverse momentum. The strategy for this 
analysis and the selection cuts to be applied were studied with a realistic 
and detailed simulation of the detector geometry and response, including the 
main background sources f3J |^] . 

The expected performance for central Pb-Pb (b < 3.5 fm) at ^/s NN = 
5.5 TeV and pp collisions at y/s = 14 TeV, as obtained using the input pro- 
duction yields Np h _ Ph = 115 and A 7 "™ = 0.16 (see Ref. [IH!), is summarized in 
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Figure 3: Double differential cross section per nucleon-nucleon collision for 
D° production as a function of p t , as it can be measured with 10 7 central 
Pb-Pb events (left) and 10 9 pp minimum-bias events (right). Statistical 
(inner bars) and p t -dependent systematic errors (outer bars) are shown. A 
normalization error of 11% for Pb-Pb and 5% for pp is not shown. 



Fig. El The accessible p t range is 1-14 GeV/c for Pb-Pb and 0.5-14 GeV/c 
for pp. In both cases the statistical error (corresponding to 1 month of data- 
taking for Pb-Pb and to 9 months for pp) is better than 15-20% and the 
systematic error (acceptance and efficiency corrections, subtraction of the 
feed-down from B — > D° + X decays, cross section normalization, centrality 
selection for Pb-Pb) is better than 20%. More details are given in Ref. |14j . 



5 Results: Raa an d Rq/h 

The nuclear modification factor for D° mesons is reported in Fig. HJ Nuclear 
shadowing, parton intrinsic transverse momentum broadening and energy 
loss are included. The dead cone effect is not included in the left-hand 
panel and included in right-hand panel. Different values of the transport 
coefficient are used for illustration; we remind that the value expected on 
the basis of the pion quenching observed at RHIC is q = 4 GeV 2 /fm. The 
reported statistical (bars) and systematic (shaded area) errors are obtained 
combining the previously-mentioned errors in Pb-Pb and in pp collisions and 
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Figure 4: Nuclear modification factor for D° mesons with shadowing, intrin- 
sic k t broadening and parton energy loss. Left panel: without dead cone 
correction; right panel: with dead cone correction. Errors corresponding to 
the curve for q = are shown: bars = statistical, shaded area = systematic. 
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Figure 5: Ratio of the nuclear modification factors for D° mesons and 
for charged hadrons. Left panel: without dead cone correction; right 
panel: with dead cone correction. Errors corresponding to the curve for 
q = 0.05 GeV 2 /fm are shown: bars = statistical, shaded area = systematic. 
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considering that the contributions due to cross section normalization, feed- 
down from beauty decays and, partially, acceptance/efficiency corrections 
will cancel out in the ratio. An uncertainty of about 5% introduced in the 
extrapolation of the pp results from 14 TeV to 5.5 TeV by pQCD is also 
accounted for (see Ref. [Tl]). 

The effect of shadowing, clearly visible for q = (no energy loss) as a 
suppression of Raa, is limited to p t < 6-7 GeV/c. Above this region only 
(possible) parton energy loss is expected to affect the nuclear modification 
factor of D mesons. 

For q = 4 GeV 2 /fm and no dead cone, we find Raa reduced, with respect 
to 1, by a factor about 3 and slightly increasing with p t , from 0.3 at 6 GeV/c 
to 0.4 at 14 GeV/c. Even for a transport coefficient lower by a factor 4, 
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q — 1 GeV 2 /fm, _Raa is significantly reduced (0.5-0.6). When the dead cone 
effect is taken into account, the -Raa reduction due to quenching is found to 
be lower by about a factor 1.5-2.5, depending on q and p t . 

We point out that the estimated systematic uncertainty of about 18% may 
prevent from discriminating between a scenario with moderate quenching 
and negligible dead cone effect (e.g. q = 1 GeV 2 /fm in the left-hand panel of 
Fig. HJ) and a scenario with large quenching but also strong dead cone effect 
(e.g. q = 4 GeV 2 /fm in the right-hand panel). 

The comparison of the quenching of charm-quark-originated mesons and 
massless-parton-originated hadrons will be the best suited tool to disen- 
tangle the relative importance of energy loss and dead cone effects. The 
D / 'charged hadrons (D/h) ratio, defined as in (J2J), is presented in Fig. for 
the range 5 < p t < 14 GeV/c. We used R\ A calculated as previously de- 
scribed and i2°A) without and with dead cone, as reported in Fig. |U Being 
essentially a double ratio Pb-Pb/Pb-Pb x pp/pp, this observable is par- 
ticularly sensitive, as many systematic uncertainties cancel out (centrality 
selection and, partially, acceptance/efficiency corrections and energy extrap- 
olation by pQCD). The residual systematic error is estimated to be of about 
10-11%. 

We find that, if the dead cone correction for c quarks is not included, 
Rn/h is essentially 1 in the considered p t range, independently of the value 
of the transport coefficient, i.e. of the magnitude of the energy loss effect. 
When the dead cone is taken into account, R^/h is enhanced of a factor 
strongly dependent on the transport coefficient of the medium: e.g. 2-2.5 for 
q = 4 GeV 2 /fm and 1.5 for q = 1 GeV 2 /fm. The enhancement is decreasing 
with p t , as expected (the c quark mass becomes negligible). 

The Rj)/h ratio is, therefore, found to be enhanced, with respect to 1, 
only by the dead cone and, consequently, it appears as a very clean tool to 
investigate and quantify this effect. 

Since hadrons come mainly from gluons while D mesons come from (c) 
quarks, the D/h ratio should, in principle, be enhanced also in absence of 
dead cone effect, as a consequence of the larger energy loss of gluons with 
respect to quarks. Such enhancement is essentially not observed in the 
obtained Rn/h because it is 'compensated' by the harder fragmentation of 
charm quarks with respect to light quarks and, particularly, gluons. With 
z the typical momentum fraction taken by the hadron in the fragmenta- 
tion, ph adron = zp\ wton , and AE the average energy loss for the parton, 

/ parton \/ parton a 7-1 1 

{Pt ) = Pt ~ AE, we have 

^hadrony = ^hadron _ g ^ ^ 

meaning that the energy loss observed in the nuclear modification factor is, 
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indeed, zAE. We have, thus, to compare z c ^d AE C to % U on^hadron AE g \ UOIl . 
With 2 g iuon^hadron ~ 0.4, z c ^b ps 0.8 for p t ' > 5 GeV/c and A£ c = 
A£? g i uon /2.25 (without dead cone), we obtain 

^c— »D AE C ~ 0.9 ^gluon— >hadron A£<gi uon . (9) 

This simple estimate confirms that the quenching for D mesons is almost the 
same as for (non-charm) hadrons, if the dead cone effect is not considered. 

The errors reported in Fig. show that ALICE is expected to have 
good capabilities for the study of Rr>/h- m the range 5 < p t < 10 GeV/c 
the enhancement due to the dead cone is an effect of more than 3 o for 
q > 1 GeV 2 /fm. The comparison of the values for the transport coefficient 
extracted from the nuclear modification factor of charged hadrons and, in- 
dependently, from the D/ charged hadrons ratio shall provide an important 
test for the coherence of our understanding of the energy loss of hard probes 
propagating in the dense QCD medium formed in Pb-Pb collisions at the 
LHC. 
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